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Abstract 

In this contribution to the study of one dimensional point potentials, we prove that if we take the 
limit 5 —>■ 0 on a potential of the type vo5{y) + 2v\5'{y) + 'WoS{y — q) + 2wi5'{y — q), we obtain a 
new point potential of the type uo5{y) + 2ui5'{y), when uq and ui are related to vo, vi, wo and wi 
by a law having the structure of a group. This is the Borel subgroup of SL 2 (R). We also obtain the 
non-abelian addition law from the scattering data. The spectra of the Hamiltonian in the decoupling 
cases emerging in the study are also described in full detail. It is shown that for the ui = ±1, wi = ±1 
values of the S' couplings the singular Kurasov matrices become equivalent to Dirichlet at one side of 
the point interaction and Robin boundary conditions at the other side. 


1 Introduction 

One dimensional models with point interactions [T] have recently received much attention. They serve 
to modeling several kinds of extra thin structures nia or point defects in materials, so that effects like 
tunneling are easily studied. They are also interesting in the study of heterostructures, where they may 
appear in connection to an abrupt effective mass change [4j. The general study of point interactions of 
the free Hamiltonian Hq = — ^ is mainly due to Kurasov [51 [6] and it is based on the construction 

of self adjoint extensions of symmetric operators with identical deficiency indices. More recently Asorey, 
Mufioz-Castaheda and coworkers reformulated the theory of self adjoint extensions of symmetric operators 
over bounded domains in terms of meaningful quantities from a quantum field theoretical point of view, see 
Refs, [a[HE! and references therein. This new approach allows a rigorous study of the theory of quantum 
fields over bounded domains and the quantum boundary effects. Here, it is well known the existence of a 
four parameter family of self adjoint extensions. Some of these extensions are customarily associated to 
Hamiltonians of the type Hq plus an interaction of type a5{x — d) bS'(x — d), where S' is the derivative 
of the Dirac delta while a, b and d are fixed real numbers. This type of perturbation has interest both 
mathematical and physical in Quantum Mechanics and has been largely discussed [n uniini na da nn 

da da HI!- 

Although it is totally clear which self-adjoint extension corresponds to the ID-Dirac delta S{x — d), 
there is no consensus on which one should be assigned to its derivative 6' {x — d) [la dll dB da dH da dll • 
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Self adjoint extensions of Hq are characterized by matching conditions at x = d. The set of self adjoint 
extensions of Hq is given by the unitary group U(2) and therefore is a 4-parameter family of operators 
(see Refs. [H H]). In previous works, we have characterized perturbations of type a8{x — d) + bS'{x — d) 
by suitable matching conditions as a two parameter family of self adjoint extensions (see Refs. miin]). 
Combining this family of point interactions with some other potentials (or even with mass jumps), we have 
obtained some physical features such that transmission and reflection coefficients, bound and antibound 
states and resonance poles pT] . 

In other physical context, scalar QFT on a line, point potentials are useful to model impurities and/or 
providing external singular backgrounds where the bosons move, see e.g. |22j . The spectra of Hamiltonians 
with 5 and 5' point interactions provide one-particle states in scalar (1 -|- I)-dimensional QFT systems, see 
Refs [a 019]. In particular, configurations of two pure delta potentials added to the free Schrodinger 
Hamiltonian have found applications to describe scalar held fluctuations on external backgrounds, see 
e.g. [13121123, as the corresponding scattering waves. The same conhguration of delta interactions is 
addressed in Reference [TH] as a perturbation of the Salpeter Hamiltonian. Moreover, according to the 
idea proposed by several authors, delta point interactions allow to implement some boundary conditions 
compatible with an scalar QFT dehned on an interval, see |22| and References quoted therein. Other than 
these generalized Dirichlet boundary conditions were discussed in [20] where it is shown that the use of 
5-5' potentials provides a much larger set of admissible boundary conditions. The TGTG-formula was 
subsequently applied to compute the corresponding quantum vacuum energies between two plane parallel 
plates represented by a 5-5' potential in arbitrary space-time dimension. 

To be more specihc, the coupling a to the <5 potential mathematically describes the plasma frequencies 
in Barton’s hydrodynamical model |26| characterizing the electromagnetic properties of the conducting 
(inhnitely thin) plates. The physical meaning of the b coupling to the 5' interaction in the context of Casimir 
physics has been unveiled only recently in m- it describes the response of the orthogonal polarizability 
of a monoatomically thin plate to the electromagnetic field. 

Since the zero range potentials mimick the plates in a Casimir effect setup it is interesting to consider 
a three-plate configuration and investigate the Casimir forces by allowing to move freely the plate in 
the middle, see e.g. [H] for an introduction to Casimir Pistons, and im E3 El 1321 for recent results. 
In particular, it is meaningful to displace the central plate towards one of the other two placed in the 
boundary; thus, we find the main physical motivation to study the particular situation where two 5-5' 
interactions are superimposed. 

This is the first problem to be analyzed in this work. The outcome is surprising: a non-abelian addition 
law emerges which corresponds to the Borel subgroup of the 5'L2(K) group. The other focus of interest is 
the case where the ^'-couplings are exceptional, i.e., those couplings for which the transmission coefficients 
are zero and the plates become completely opaque: the left-right decoupling limit. It happens that the 
distinguished Dirichlet/Robin boundary conditions are implemented in this case by the <5-i5'-potentials but 
the superposition law just mentioned becomes singular. We shall fully discuss this much more awkward 
regime in the second part of the paper. 

Here, we shall combine two point potentials of the type ad5(x — d) 2b45'{x — d) with d > 0, and 
ao5{x) 2bo5'{x). In fact, the idea of using double delta potentials (6^ = bo = 0) has a long tradition. 
For instance, in condensed matter physics, in the BCS model [53, in Bose Einstein condensates |M] . 
or potentials with double pole resonances |3S| arise in the description of some unstable states. Finally, 
regarding physical contexts where 5-5' interactions play a role, we mention that point supported potentials 
can be used to model gap-like impurities in graphene layers and nano-ribbons, see Ref. |36| . In this 
physical situations graphene surface plasmons suffer total reflection after collision with a point-supported 
impurity/edge. 

In the present article and after the introduction of our Hamiltonian in Section, 2, we study bound states, 
scattering coefficients and resonances in the so called regular cases, which are those with b^ ^ ^ bo, 

which is done in Section 3. Next in Section 4, we consider the limit d —>■ 0. In this limit both interactions 
coincide. The result is again an interaction of the type a5{x) -\- b5'(x) and it is quite interesting to 
evaluate the composition law for the coefficients, which is not linear as one may expect naively. Indeed 
this composition law establishes the two-dimensional space of couplings (a, b) ^ as the Borel subgroup 
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of 5'L2(M), a quite unexpected result. Section 5 is devoted to study the interesting cases of left-right 

fc2 

decoupling limit, i.e., those with 6^ = and &o = The paper is closed with some concluding 

remarks. 


2 The Hamiltonian 

Let us consider a one dimensional free Hamiltonian Hq with a potential of the type a5{x — d) + bS'(x — d). 
Its Schrodinger equation reads, 

- - + adix — d)il){x) + bS'ix — d)4’{x) = Etjjix). (2.1) 

2m dx‘‘ 

In order to work with dimensionless quantities, let us introduce new variables and parameters 

h h 2a mb 2E / \ \ 

x = — y, d= — q, wo = wi = —, e = — (p{y) = tp{x), (2.2) 

me me he m mc‘^ 


such that (2.1) becomes 


dy‘ 


‘P{y) + WoS{y - q)(p{y) + 2wi6'{y - q)(p{y) = e (p{y). 


(2.3) 


From now on, we will consider this version of the Schrodinger equation instead of 

The point potential we are interested in, woS{y — q) + 2wi5'[y — q), is usually defined via the theory 
of self-adjoint extensions of symmetric operators of equal deficiency indices so that the total Hamil¬ 

tonian iJ = TLq + woS(y — q) + 2wi6'{y — q) is self adjoint. The crucial point is finding the domain of 
wave functions f{y) that makes E[q self adjoint over the domain K/jg} and characterizes the potential 
woS{y — q) + 2wi5'{y — q). As these functions and their derivatives should have a discontinuity at y = q, 


we have to define the products of the form d{y — q)ip{y) and S'{y — q)ip{y) in (2.3). These can be done in 
several ways [nmaiiT], but we choose the following: 




(2.4) 


ny - yMy) = _,) _ _,, , (2.5, 


where f{q'^) and f{q ) are the right and left limits of the function f{y) as y —>■ q, respectively. The 


Schrodinger equation (2.3) should be viewed as a relation between distributions. 


In order to obtain a self-adjoint determination of the Hamiltonian El = E[q + wod{y — q) + 2wiS'{y — q), 
we have to find a self adjoint extension of EIq. In order to do it, we have to find a domain on which 
this extension acts. This domain is given by a space of square integrable functions satisfying certain 
assumptions including matching conditions at the point q that affects to the value of wave functions and 
their derivatives at q DEI- In particular, this implies that both wave functions and derivatives cannot be 
continuous at q, so that equations (2.4) and (2.5) make sense. 

The functions in the domain of the Hamiltonian El are functions in the Sobolev spac^ IT|(M \ {q}) 

^This is the space of absolutely continuous functions f(y) with absolutely continuous derivative f'{y)y both having arbitrary 
discontinuities at q, such that the Lebesgue integral given by 


r 


{\f(y)\^ + \f"iy)\^}dy 


converges. 
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such that at q satisfy the following matching 


condition^ 


V ‘^'(9^) 


1 + Wl 

1 — Wl 

Wq 

1 — wf 


0 

1 — Wi 
1 + Wl 


V5(<? ) \ 
^'(9-) ) ■ 


( 2 . 6 ) 


These results could be, in principle, extended to interactions of the type o,iS{y — 9i) + ‘2.biS'{y — qi), 
where N could be either finite or infinite [S]. In the present paper, we assume that TV = 2 as above, so 
that the total Hamiltonian takes the form: 


H = Ho + V + W , 


Ho = - 


dy"^ ’ 


(2.7) 


with 

V = vo5{y) + 2vi5'{y), W = wo5{y - q) + 2wi5{y - q), (2.8) 

where Vq = ^, Vi = , and A, B are respectively the <5 and 5' couplings of the pair placed at the origin 

in the original variables. While V is supported at the origin ?/ = 0, W is supported at a point q that we 
assume positive, g > 0. The potential V + W \s physically relevant as is related to the Casimir effect [5D]. 
The corresponding dimensionless Schrodinger equation is 


d'^‘p{y) 

dy'^ 


bo^(y) + ‘^viS'iy) + wo5{y - q) + 2wiS'{y - g)] (p{y) 


e<p(2/)- 


(2.9) 


One of the motivations of the present work was the study of the effect resulting of taking the limit in (2.9) 
as g —>■ 0, i.e., the effect of the superposition of the two point potentials V and W at the same point. We 
shall see that, as a result of the limit, we obtain a point potential of the type uoS{y) + 2ui5'{y), where uq 
and Ui are not the sums Vq + Wq and vi+Wi, but instead another kind of law, which has a group structure 
as the Borel subgroup of S'L 2 (K). Some further results will be discussed. 


3 Matching conditions and scattering coefficients 



woS(y-{ 

r) + 2u.'i S'{y - 

region 1 

region 2 

region 3 

0 


9 

voHy)- 

f 2i', S'{y) 



Figure 1: Schematic representation of the physical problem under study, given by ( |2.9[ ), before considering the 
limit g —>■ 0: a potential with a couple of singular point interactions described by a linear combination of Dirac 
delta distributions and their derivatives. 


In order to analyze the solution of (2.9), let us split the real line into three regions, as is shown in 


Figure [2 region 1 is the half line with ?/ < 0, region 2 is the segment 0 < y < q and region 3 is the half line 
y > q. Supported on y = 0 and on y = g, we have the point potentials V and W, as in (2.8), respectively. 


In these regions the values of the solution and their derivatives are those of the free Schrodinger equation, 
given by: 


^j{y) = Aje-^'^y + Bje^'^y, ^'Jy) =-ik{A^ - B^ , j = I,2,3, = e > 0. 


(3.1) 


^This is true for wi 7 ^ ±1, while for wi = dil, we have to define the matching conditions in another way [^. We shall 
concentrate now in the regular cases, but the decoupling limit will be also considered later on in Section 
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Since the point potential V is defined by matching conditions like those in (2.6), we have that 


A 2 + 

—ik{A2 — B 2 ) 


Defining the following matrix: 


= M„ 


Ai + Bi 

—ik{Ai — Bi) 


M„ 


^ 1 + ^1 
1-Vi 

Vo 


0 

1-Vi 


\ 1-vl 1 + 




K = 


1 1 
—ik ik 


equation (3.2) becomes: 


A 2 

B 2 


= K-^ M„ K 


Ax 

Bi 


(3.2) 


(3.3) 


(3.4) 


This expression (3.4) gives the matching conditions at the point y = 0. At the point y = q the same 

(3.5) 


procedure works after appropriate translation 

= Q-^K-^M^KQ 


and we finally obtain: 


A 3 

B 3 


A 3 

B 3 


= Q-^ K-^ K Q My K 


A 2 

B 2 


-1 


Ai 

Bi 


= T„ 


Ai 

Bi 


(3.6) 


where the definition of Tq is obvious, while Q and My; are respectively: 


Q = 


„ — iqk 


oiqk 


1 + 

M - I 

IVly, — I 


0 

1 — Wi 

1 — wf 1 + 


(3.7) 


The Tq-matrix in (3.6) relates the asymptotic behaviour of the two linearly independent dost scattering 
solutions at a: << 0 with their counterparts at a: >> 0, see e.g. m- The Tq-matrix elements satisfy the 
identities: 

J _ rp']_2rp']_2 _ 1 ^11 _ rTi22 rTil2 _ rTi21 

~ -^q ^ q — q ’ q — q 

such that Tq is in general an element of the group SL 2 (C). From the T^-matrix one obtains the scattering 
matrix Sq by an standard procedure. Reshuffling the linear system (3.6) in the form 


Ai 

B 3 


= Sy 


A 3 

Bi 


one easily checks that 

1 


S„ = 


Til [ T 


1 

q 


detTo 


5t = J- 

9 7^11 

q 


1 


rp2l 

-fl^ AetTq 


SlSn = f 


-'q'-’q 


0 1 


(3.8) 


(3.9) 


i.e., the scattering ^g-matrix arises as a 2 x 2 unitary matrix: = SqS}^ = 1. The transmission and 

reflection coefficients provide the usual form of writing the S'^-matrix elements: 


Sq = 


t{k) rii{k) 
rhik) t{k) 


(3.10) 


t{k) is the amplitude of the transmitted waves coming from the far left or from the far right. Time-reversal 
invariant potentials give rise to identical transmission coefficients: tfi{k]q) = t]^{k\q) = t{k\q). However 
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the reflection amplitudes are different for incoming waves either from the far right or the far left such that 
‘>'R{k',q) 7 ^ TL{k\q) because the potential is not parity invariant. Comparison of equations (3.91 and (3.10) 
shows that: 

1 ri2(fc) T2i(fc) 

= ^11/, N . rR{k\q) = - , rL{k,q) = ■ (3.11) 

r“(fc) 

Transmission and reflection coefficients were calculated in |20j when the two 5-5' point interactions were 
symetrically located with respect to the origin. For the arrangement chosen in this paper we find from 
formulas (13.Ill) : 


rL{k-,q) = - 


rRik;q) = 


t{k;q) = 


e (2k (uj + l) + ivo) {-ikwi — iwo) + {4kvi — ivo) {2k {wf + l) — iwo) 

W) ^ 

e 2 iqk ^2k {vf + l) — ivo) (4kwi + iwo) + {4kvi + ivo) {2k {wf + l) + iwo) 


A(fc) 


4fc^(l — Ui)(l — wf) 

W) 


A(fc) = + 4:ikvi){wo — Aikwi) + {2kvf + 2fc + ivo){2kwf + 2fc + iwo), 

where A(fc) is a function of k and the other parameters of the problem. 

3.1 Bound/antibound states and resonances 

Complex zeroes of 


(3.12) 

(3.13) 

(3.14) 

(3.15) 


T^\k) = 


A(fc) 


(3.16) 


4fc2(l — r'i)(l — wf) ’ 

which are poles of the Sg-matrix in the /c-complex plane, give rise to bound or antibound states if are 
located on the purely imaginary, respectively positive or negative half-axis. Complex zeroes coming in 
pairs having opposite real part and identical imaginary part correspond to resonances. 

The relevant physical information is encoded in the analysis of complex zeroes of T^^(k), or equivalently 
of A(fc). In order to simplify this analysis, it is useful to make the following changes in the parameters in 

(P^: 

qwo Vi 


2kq = z, a 
Then, A(fc) = 0 is equivalent to 


qvo 

l + v^ 


Wi 




2 ’ 


1 


1 -I- uif 


z ia \ 

a -\- 2ivz J 


z ir \ 
T — 2iwz 




Z = Zr iZi, Zr,Zi G 


(3.17) 


(3.18) 


Equation (3.18) can be considered as a complex or two-dimensional generalization of the so-called gener¬ 
alized Lambert equation (see |38L I39j and references quoted therein). From this complex Lambert equation 
we can obtain two real equations, corresponding to the real and imaginary parts of ( |3.18 ): 

[4:vwz‘^ — {2vZi — a){2wzi -\- r)] cosz^ ~ 2Zr{TV -\- AvwZi — aw) sinz^ = [{zi + a){zi -I- r) — z^] (3.19) 

2Zr{TV -\- AvwZi — aw) cos Zr -h [4vwz^ — (2vZi — a)(2wzi -I- t)] sinz^ = —Zr(2zi -|- cr -I- T)e^'. (3.20) 

These expressions for real and imaginary parts of A(fc) = 0 are rather complicated. A simpler compatibility 
condition is given by 

2z, + {^VZ, - cr)2) (dw^Z^ -t (2wZi -k t)^) 

which is, indeed, a generalization of the Lambert equation that includes two real variables. The solutions 
we are looking for must satisfy (3.21), although (3.19)-(3.20l are more restrictive. First of all, we can 
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easily check that a simple solution of the system (3.19)-(3.201 is = 0, Zi = 0, but this corresponds to 
k = 0, which is a pole of T^^{k). In addition, it is obvious that the system is symmetric on the variable 


Zr, although has no symmetry on Zi. The typical behavior of the solutions of (3.19)-(3.21) is shown in 
Figure 




Re(z) 


— Curve (3.19) 

— Curve (3.20) 

— Curve (3.21) 


Figure 2: On the left, for some values of the parameters {vo = —2, vi = 4, wq = —1, wi = 3 and q = 1/2), the 
typical behavior of the curves (3.191-( 3.211 is plotted. The intersections between the blue and orange curves (the 
values of z such that the real and imaginary parts of A(z) are simultaneously zero) are seen as discrete points, 
placed on the compatibility condition curve (green). Remark that there is an inhnity number of such solutions on 
the lower part of the complex plane, appearing in pairs z = ±2r,n — izi,n, Zr,n, Zi^„ > 0, n = 1, 2,..., corresponding 
accordingly to resonances in the spectrum; in addition, there is also one imaginary solution with positive imaginary 
part zo = izifi. On the right, a zoom of the same plot near the origin is shown. One may observe a bound state in 
the upper intersection. The lower intersection, however, happens exactly at the origin and it is not a meaningful 
physical solution. 


4 The addition law at the q = 0 limit 

Proposition 1. The Tq matrix for two-pairs of 6-S' interactions displaced from each other a distance q 


is determined in formula (3.6). In the limit q ^ 0 Tq becomes: 

= ■ M„K ‘ 


^3 

^3 


^1 

Bi 


= R-^ M„K 


Ai 


(4.1) 


That is to say, the superposition of these two S-S' is equivalent to a single 5-5' where the matching conditions 
are characterized by the Rurasov matrix 


l + Ml 


M — ^““1 

- Un 


0 


Uq 1 — Ui 

1 — U^ 1 + ltl 


Ul+Wi Uo(l - Wl)^ + Uio(l + Ui)2 

where ui = -—^-, uq = -——^- — -. (4.2) 


1 + ViWi 


(1 + ui'u;i)2 


Proof. In the q = 0 limit, the transfer matrix is lim,j_>o Tq = Tq, while Q obviously becomes the identity 
matrix. Since and My do not depend on q, the q = 0 limit of (3.6) reads: 


As 

Bs 


= R~^ My, My R 


Ai 

Bi 


= Tq 


Ai 

Bi 


(4.3) 
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One may expect that the product of two Kurasov matrices My is another Kurasov matrix. This 
assumption means that 


My, My 


1 + 

1 — Wi 

Wq 

1 — w\ 


0 

1 — Wi 
1 + 


0 

1 — Vi 
VQ 1 - 

\ l-vf 1 + Ui 


/1^ Q 

1 — Ui 

Uo 1 — Ui 

\ 1 — ul 1 + Ul 


= My. 


(4.4) 


A lengthy but straightforward calculation shows that this is the case with the “composite”couplings: 


vi+wi _ ?;o(l - 

l + uiici’ (l + ui'u;i)2 


(4.5) 


Thus, the collapsed two S-6' pairs are tantamount to a single point potential of the form UQ6{y) +2uiS'{y) 
defined by the transfer matrix Tg = K~^MyK, q. e. d. 

One would have expected another result such as an additive rule of the type ug = fg + wq and ui = 
vi + wi, due to the form of point potentials that converge. In fact for zero S' couplings we would obtain: 




and the fusion of two S point potentials is a pure abelian process. Nonetheless, the addition law (4.5) is 


quite interesting. First of all, the first equation in (4.5) resembles the addition law of velocities in special 
relativity. However, the second one looks more complicated. In any case, one natural question we may 


pose is if the composition law (4.5) has a structure of a group and, if this is the case, which group this can 
be. 

Checking the group structure is quite simple. The product of two matrices My and My, as in (4.4) 


gives another matrix with identical structure. This gives the product law. The associativity comes from 
the associativity of the product of matrices. The identity of the group is the 2x2 identity matrix I, which 
corresponds to take both parameters equal to zero in (2.6), Mg = I. The inverse of M„ is a matrix M, 


so 


that My Ms = I. The calculation is straightforward. For My as defined in (3.2) one finds 


Ms = 


( ^ ~ '^1 
1 + r"! 

-Ug 


0 

1 + fl 


(4.6) 


\ 1 — vf \ — Vi 


which shows that taking the inverse in our group of matrices is equivalent to the transformation (ug, Ui) —)■ 
(—ug, —rii). We may denote the inverse of My as M_„. Thus, the structure of group in our set of matrices 
has been confirmed. Furthermore, there are two other remarkable properties of matrices in this group: 

(i) Non commutativity: 


0 

MyMy, MyyMy = | 4('Cg'lCl - tCgUl) 


(l-n2)(l- 


wl) 


(ii) The trace of M„ is independent of vg: 


tr My = 2 


1 + ^1 
1-vl 


(4.7) 


(4.8) 


This trace has two singular points Vi = ±1, emerging from Kurasov’s analysis, and corresponding to the 
decoupling limits which are not included in the type of matrices that we are considering in the present 
section (for an analysis of the decoupling or completely opaque limits see Section]^. It is positive on the 
interval —1 < ui < 1, reaching the minimum value 2 for vi = 0. On the half lines —oo < rii < —1 and 
1 < Ul < 00 the trace is negative and its modulus is bigger than 2. Its limits when vi —> ±oo are equal to 
—2. A plot of the function (4.8) can be seen on Figure]^ 





























Figure 3: 


Plot of the trace of a generic matching matrix M„, 


given by equation (4.8 1 . 


4.1 The Borel subgroup of SL 2 (M.) and its Lie algebra 


Proposition 2. 

group law: 


The Kurasov Matrices form a group where the matrix product My ■ My, = My is the 


vi + wi _ t;o(l - wi)^ + wo{l + ui)^ 

l + UlWi’ (l+UiUii)2 


(4.9) 


The group of Kurasov Matrices is the Borel subgroup of SL 2 {M.) 

Proof. All the previous information shows that the set of My matrices is the subgroup of SL 2 {M.) of all 
the lower triangular matrices such that |tr(( 7 )| > 2. This is the Borel subgroup of SL 2 {M.). 

The Lie algebra of the Borel subgroup of SL 2 {^) is the minimal parabolic subalgebra of 5(2 (M). The 
Cartan basis for sl 2 (K) is given by the following 2x2 matrices: 



Their commutation relations are 




[e, f] = h, [h, f] = - 2 / , [h, e] = 2 e. 


(4.10) 


(4.11) 


The Borel subalgebra of a given Lie algebra g with Cartan subalgebra t) is the direct sum of the Cartan 
subalgebra with a given ordering and the weight spaces of g with negative weight (it is equivalent to take 
positive weights; this would just change upper diagonal matrices by lower diagonal matrices. In any case, 
both algebras, or equivalently both groups, are isomorphic). Therefore in the case of sl 2 (K) the Borel 
subalgebra is 

9KM = ihJ)- (4.12) 

The elements of the two disconnected components of the group are obtained by exponentiation: 


^ah+Pf _ 


P 


gQ 

sinha 

a 



—e 


ah-\-l3f _ 


-P 


—e 

sinh a 
a 


(4.13) 


For the component of the elements connected with the identity, the couplings are easily obtained in terms 
of the group parameters: 


^ah+Pf _ 


= My 


^ P , OL a 

vq = 2 — tanh—, vi = tanh— 
a 2 2 


(4.14) 
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On the other hand, for the elements in the other connected component, we have: 

_e°'h+Pf — vq = 2 — coth—, vi = coth— . 

a 2 2 

The composition law of the group can be expressed in terms of these exponentials. In fact, 

^aih+/3if ^a2h+l32f _ ^ah+/3f 


with 


a = ai + a 2 : 


/? = 


(ai + 02 ) (e“^Q;2/3i sinh (ai) + e “^01/32 sinh ( 0 : 2 )) 


aiQ ;2 sinh(Q!i + 0 : 2 ) 
This concludes the discussion on the group structure. Q. E. D 


(4.15) 


(4.16) 

(4.17) 


4.2 Reflection and transmission coefflcients due to a single S-S' interaction at 
the origin 

Proposition 3. The q = 0 limit of the scattering coefficients of two a priori separated pairs of 5-5' 
interactions exactly coincide with the scattering coefficients produced hy a single 5-5' pair with couplings 
determined from the Kurasov matrix. In terms of the scattering matrix we can write: 


lim = S{My ■ M„) = S{Mu). 


(4.18) 


Hence the scattering matrices produced by two superimposed pairs of 6-5' interactions give rise to a repre¬ 
sentation of the non-abelian composition law of Kurasov matrices. 

Proof. The linear map between asymptotic scattering dost solutions respectively in the far right and 
the far left. 



produced by a single point interaction V{x) = uq5(x) -\-2ui6'{x) is provided by a similarity transformation 
of the Kurasov matrix : 


Z = K-^M^K , 


We obtain for the .Z-matrix elements 


1 + Ul 
1 — Ul 

Uo 

l-uj 


1 - Ml I 

1 + Ml / 


Z^Hk) 


2fc(l + ul) + iup _ ^22 iu\ 
2k{l-ul) 


Z^\k) = 


4kui + iuQ 
2fc(l — Ui) 


= Z 


21 


(4.19) 


where the bar stands for complex conjugation. From this information one reads the scattering coefficients 


t{k) 


Z^'^fk) 4A:rti + luo 

Z^i(fc) 2k{ui + 1) + iuQ ’ 
1 _ 2k{ul - 1) 

Z^^fk) 2k(u\1)-\-iuf)'' 


, , Z^^fk) 4:kui — iuo 
^ Z^^{k) 2k(u\-\-1)-\-iuQ^ 


(4.20) 


(4.21) 


in perfect agreement with the results in m- 

Finally, we compare these scattering coefficients with those obtained from the Tg-matrix that describe 
the map between asymptotic dost solutions when the distance q from one 5-5' point interaction to the other 
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tends to zero. The outcome is: 


(4.22) 

(4.23) 

(4.24) 

to check that the addition law 

vi+wi ?;o(l-Wi)^+u>o(l + Ui)^ 

1 I ’ /'I I ^2 ’ 

1 + ViWi (1 +'Ci'lCi)"^ 

also works for the scattering coefficients, q.e.d. 


( . _ TjJ^(A:) _ vi ({vi + 2) two + 4ifc) + AikviWi {vi + wi) + 4ikwi + vq (wi — 1) ^ + icq 

^ ’ Tg^{k) —2ikwi{{v\ + \)wi + Avi) — 2ik{v‘l-\-l) + {vi + 1 )'^Wq + Vq{wi — 

, , _ T§^{k) _ vi (— (t;i + 2) Wo + 4i/c) + 4:ikviWi {vi + wi) + Aikwi — vq {wi — 1) ^ — wq 

’ T^^ik) -2i/i:wi (( 1 ;^ + 1) wi + 4r;i) - 2z/c (t;^ + 1) + (ui + 1) 2^0 + uo - 1) ^ 

t(k-O) = ^ = _ 2fc {vf - 1 ) {wj - 1 ) _ 

T}^{k) 2kwi[{vl+ l)wi +'ivi)+ 2k{vl+ 1) + i(vi + l)‘^WQ + ivQ(wi - l)'^' 


5 The left-right decoupling values of the 5' couplings 


We have seen that there are singularities of the Kurasov matrix at the critical points, ui = ±1 and 
wi = ±1. One may expect that these critical cases where the transmission coefficients are zero do not 
contribute to the group structure that we have previously analyzed. The completely opaque potentials as 


given in equations (2.8) have clearly one of the following four forms: 


V = voS(x) ± 2S'(x) and/or W = wo6{x — g) ± 23'(x — q). 


(5.1) 


In order to define these potentials, we obviously cannot use matching conditions of the form (2.6) that are 
singular. Instead, we shall impose [S]: 


1. For vi = 1: 


2. For wi = 1: 


)=0, (p'(0+) = ^(/2(0+). 


(p{q ) =0, (fi'iq^) = ^ifiq^). 


(5.2) 


(5.3) 


Thus, we set Dirichlet boundary conditions on the left of the two points a: = 0 and x = q, whereas 
Robin boundary conditions are chosen on the right hand side of a: = 0 and x = q, see |2()j . In the 
two remaining cases the matching conditions are the mirror images of the previous ones. 


3. For vi = —1: 


(^(0+)=0, (p'{0 ) = ~ip{0 ). 


(5.4) 


4. For wi = —1: 


Wo 


ip{q+) = 0, ip'{q ) = —-ifiq ) . 


(5.5) 


We remain, however, in the physical situation depicted in Figure were the plane waves and their deriva¬ 
tives are written as in (3.1). There are eight possible configurations involving at least one decoupling 


configuration of the couplings: either the two S-S' interactions build opaque walls both at a; = 0 and x = q, 
or, there is no transmission only at one point. We discuss first the cases when the two S' couplings take 
the decoupling limit. 
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5.1 Two 6' couplings in the decoupling limit 


There are four cases in which we have decoupling situations both at x = 0 and x = q. Let us consider 
them separately. 

Case 1 : vi = l,wi = 1. 

The boundary conditions are given by (5.2) and (5.3). Let us consider the situation on the interval [0, g]. 
For X = 0, vj'(0+) = ^ is written as: 


-ik(A2 - B 2 ) = —(A 2 +B 2 ) ^ - 

^ 4 Vo + Aik 


A 2 


vq — Aik 


Ak 


— = — exp —2i arctan — 

Vo 


For X = g, we write (p{q ) = 0 as 


A 2 + B 2 = 0 ^ — = . 

B 2 


(5.6) 


(5.7) 


Then, we have a transcendental equation, which can be written either in the form of a generalized Lambert 
equation: 

_ ., 7 7/-. — Aih 

(5.8) 


^2^kq ^ VQ - 

Vo + Aik ’ 


or in terms of the function arctan as 


7 4^ \ 4^ 

kq = — arctan — tan(Kg) =-. 

Vo 


(5.9) 


This transcendental equation has a countably infinite number of solutions, kn, that give the energy levels 
corresponding to this situation. In the limit g = 0, we have only one solution, k = 0. Outside the interval 
[ 0 ,g], we have the equations 


V3(0 ) = 0 => Ai + Bi = 0; 




(5.10) 


These are relations between coefficients, which do not provide of any further information, so that we shall 
not refer for similar situations which will appear for all other cases. 

Case 2 : vi = 1, Wi = —1. 


The boundary conditions are (5.2) for x = 0 and (5.5) for x = g. The condition at x = 0 has already been 
studied, providing equation (5.6), so that let us consider the new boundary condition at x = g. It comes 
from (p'(g-) =gj(g-): 


^2 


—e 


2^kq '^0 + 

Wq — Aik 


_ _^ 2 ikq I arctan —^ 


Ak 

Wo 


(5.11) 


This equation is to be compared to (5.6). The result given in terms of a generalized 




^2ikq ^ ^0 - 4»fc Wq - Aik 
Vo + Aik Wo + Aik 


(5.12) 


This transcendental equation can be written in another way by playing with the formulas for the arctan. 
Note that from (5.12) we obtain straightforwardly the following expression: 


Ak 

Ak 

arctan — 

+ arctan — 

Wo 

■^0 


tan(A:g) 


Ak{wo+Vo) 
wovo— 16 k^ 


(5.13) 


This new transcendental equation gives another set of energy eigenvalues (see Figure [fright). In the limit 
g = 0 , we have only one solution: k = 0. 
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Case 3 : Vi = —1, Wi = 1. 

The boundary conditions are (5.4) and (5.3). This means (/^(0+) 
A 2 /B 2 = —1 and = 0, which is 


0 , which implies that ^2 + i ?2 = 0 or 


^26 


— ikq 


^ 26 *'=® = 0 ^ = -e" 

^2 


^2 _ _ 


so that 


A 2 


= -1 = 


k = 


(5.14) 

(5.15) 


Again, in the limit g = 0, the only solution is fc = 0. 

Case 4 : vi = —1, wi = —1. 

This correspond to (5.4) and (5.5). The boundary condition at a; = 0 gives A 2 /B 2 = —1. The boundary 
condition at a: = g is just that given in (5.11). Then, the transcendental equation giving the energy levels 
is given by 

7/lri - ZL? h 

(5.16) 


^2^kq ^ 

Wo + 4ik 


in the form of a generalized Lambert equation, or 


tan(A:g) = — 


4k 

Wo 


(5.17) 


in the form of a transcendental equation in terms of the tangent. This transcendental equation also has a 
countably infinite number of solutions, fc„, that give the energy levels corresponding to this situation, (see 
Figure]^ left). Once more, in the limit <7 = 0, only the solution k = 0 remains. 


Wi=\/i =-1: ]/o=Wo=-2 Wi =- 1/1 =-1; \/o=wo=-2 




— Curve (3.19) 
Curve (3.20) 


Figure 4: Plots of the equation Re(A( 2 )) = Im(A( 2 )) = 0 over the complex 2 -plane when both singularities have 
couplings (as indicated on the top) in the decoupling regime. The green line represents the axis Im( 2 ) = 0. 
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5.2 Only one 5' coupling in the decoupling limit 

Let us consider now the four cases in which we have a decoupling and a regular coupling. 

Case 1 : vi = l,rci ^ ±1. 

From vi = 1, that is, from (5.2), we know that the conditions Ai = —Bi ^ 0 and (5.6) must be satisfied, 
that is 


(pi{y) (X sin(fcy). 


A, = B,. 


(5.18) 


Vo + 4ik 

Note that the system behaves as if there is an impenetrable barrier at a; = 0. The coefficients A 3 , B 3 are 
obtained from A 2 ,B 2 as in (3.5) 



/ 2 fc(l + wl) + iwo 
2 fc(l — wl) 
-2^k^ 4fcu;i - iwo 

V ^ ^ 


2*fc9 \ 

2A:(1 — wl) 
2 fc(l + wl) — zwo 

/ 



(5.19) 


2 A :(1 — wl) 2 fc(l — wf) 

All the relevant information about bound states, antibound states and resonances is obtained imposing the 
so-called purely outgoing boundary condition, which in our case is A 3 = 0. Using (5.18) and (5.19) we get 

Vq — 4ik 2A:(1 -|- wl) + iwo 


2ikq 


Vo + 4ik 4kwi + iwo 


(5.20) 


This equation is the equivalent of A(fc) = 0, where A(A:) is given in (3.15), which provides the relevant 
information in the non-decoupling case. Indeed, ( |5.20| ) is obtained making ui = 1 in A(fc) = 0. The 
analysis in this case is similar to the one carried out previously in Section M 
We are interested especially in the case where q ^ 0. Then, from (5.20) we get 


(uo + 4ik){4kwi + iwo) = (vo — 4ik)(2k(l + wl) + iwo)- 
The complex solutions of this equation are the following: 


fcn = 0, 


ki = —i 


. 4wo -I- uo(l — wi)^ 


4(1 -I- wi)^ 


(5.21) 


(5.22) 


The solution fc = 0 is not relevant, but the other produces interesting results: If 4wo -f uo(l — wi)^ < 0, ki 
is on the positive imaginary axis and it corresponds to a bound state, and if dzco +uo(l — wi)^ > 0 , fci is on 
the negative imaginary axis and it corresponds to an anti-bound state. The curves solving the equations 
Re(A) = Im(A) = 0 for this case are represented in Figure]^ left 

Case 2 : vi = —1, wi ^ ±1. 

The system behaves also as if there were an impenetrable barrier at x = 0. Now, from ui = — 1 (or 
equivalently from (5.4)), we find the condition A 2 IB 2 = —1 and consequently ( |5.19[ ), which must be also 
satisfied. Then, the purely outgoing boundary condition A 3 = 0 implies that (see Figure [fright) 


2ikq _ ‘ 2 .k{l + wl) + iwp 


e " = 


4kwi + iwo 


(5.23) 


In the limit g —> 0 we get the condition fc(l — zci)^ = 0, which has the unique solution k = 0. 
Case 3 : vi ^ ±1, wi = 1. 


In this situation we must take into account (5.3) and (3.4). The first pair of equations shows the presence 
of a kind of impenetrable barrier, now in x = 9 , and also the fact that A 2 = —Equation (3.4) can 
be rewritten in the form 



= K-^MZ^K 



/ 2 k{l + vl) — ivo 
2 fc(l — vl) 
—4kvi + ivo 
V 2 k{l-vl) 


4kvi + ivo \ 
2 fc(l — vl) 


2 k{l 


■IVo 


2 k{l-vl) / 



(5.24) 
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The purely outgoing boundary condition is in the present case is Bi = 0, so that 


g2ifcg _ 


2 fc(l + vf) + ivo 
Akvi — ivo 


(5.25) 


In the limit q ^ 0, we get fc(l + vi)^ = 0, which as the unique solution k = 0. Observe that this case is 
similar to the previous one, Case 2. 

Case 4 : vi ^ ±1, wi = —1. 

For wi = —1, the second equation in 


5.5) imposes ip'{q ) = — ^ ). Then 


M — -e -— B 2 . 

Wq — 4ik 


Taking this into (5.25), and imposing the purely outgoing boundary condition Bi = 0, we get 


g2ifeg _ _ 


Wo — 4ik 2k(l + Vi) + ivo 
Wo + 4ik 4kvi — ivo 

In the limit case g —>■ 0, we have the following complex solutions 

4vo + uio(l + viY 


ko = 0 , 


ki = —i 


4(l-ui)2 


(5.26) 


(5.27) 


(5.28) 


As in Case 1 studied before, the solution A: = 0 is not relevant, but the other one, ki, corresponds either 
to a bound state if drig + u;o(l + Ui)^ < 0 or to an anti-bound state if Avq + rco(l + ^'i)^ > 0- 


1^1=1, w^=2]Vq=Wo=-2 \/i=-1,Wi=2: v'o=Wo=-2 




— Curve (3.19) 

— Curve (3.20) 


Figure 5: Plots of the equation Re(A( 2 ) = Im(A( 2 ) = 0 over the complex 2 -plane when one singularity, as indicated 
on the top, takes the decoupling value for the coupling. The green line represents the axis Im( 2 ) = 0. 


6 Concluding remarks 

We have studied a one-dimensional two-point perturbation of the free Hamiltonian Ho = —d^/dx'^ of the 
form VoS(x) -I- 2viS'(x) + WoS(x — q) + 2wi6'{x — q) with g > 0. This type of potentials have interest in one 
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dimensional quantum field theory and on the study of the Casimir effect as shown by previous work of our 
group [2211111 HO]- We have obtained two remarkable unexpected results. 

The former refers to the limit q —>■ 0. The result is a point potential of the form uo6{x) + 2uiS'{x), 
where uq and ui are not the sums vq + wq and vi + wi respectively, but more complicated functions of these 
arguments instead. A law of composition for the coefficients of the deltas is produced with structure of 
group which coincides with the Borel subgroup of SL 2 (R). From the analysis of reflection and transmission 
coefficients, we obtain the same group law. 

The second one comes from the analysis of the decoupling cases which arises by the choices ui = ±1 
(decoupling case at a; = 0) and wi = ±1 (decoupling case at a; = g). When considering the decoupling cases 
at a; = 0 and x = q, the barrier seems to be impenetrable at these points so that no scattering states are 
produced. One finds bound or antibound states as solutions of a generalized Lambert equation on terms 
of the momentum. When the decoupling regime is produced at one point only (either x = 0 or a; = g), we 
have impenetrability at this point and scattering through the other one and resonances (as pair of poles 
of the analytic continuation of the S'-matrix in the momentum representation) can be found. However, it 
is quite remarkable to note that these poles may lie on the real axis if at least one of the 6 ' couplings is 
complex. This situation may violate a widely accepted causality condition [dOj . 

Finally, we mention that in the course of this paper we have clarified the relation between two apparently 
different ways of characterizing self-adjoint extensions of the Hq operator. The first one is based on the 
matching conditions determined from the Kurasov matrices. We have shown that they form the Borel 
subgroup of 5'L2(R)- The standard approach to deal with self-adjoint extensions of symmetric operators, 
coming back to von Neumann, is through unitary matrices, see e.g. the previous works mw by Asorey, 
Munoz-Castaheda et al. The connection from these two points of view starts from the isomorphism between 
the SL 2 (U.) and SUi_i{C) groups. In fact they are conjugate subgroups inside GL 2 (M): 


gSL2{M.)g — SUij 


a = 


1 —i 
1 i 


The matching conditions (2.6) become in the 5'f7iu(C) framework: 

= IF, 


ip{q+) - z:^'(g+) 
ip{q+) + iip'{q+) 


tp{q ) - iip'iq ) 
(f{q-) + i(p'{q-) 


where the IF,, is conjugated to the Kurasov matrix My through the action of g: 

1 / 0^1 , .,2 


Wy = g My g = 


2 ( 1 -u?) 


2(1 + vf) — ivo 4vi — ivg 

4ui + ivo 2(1 + Vi) + ivo 


( 6 . 1 ) 


( 6 . 2 ) 


In a parallel reshufling of the linear system (6.1) to that performed to dehne the Sg scattering matrix from 
the Tq matrix and passing from the (3.6) equation to (3.8) and (3.9) we rewrite (6.1) in the form 


( )- w'{.a ) \ ^TT ( ‘/5(g+) - \ 

V + V(9+) ) V + ) ’ 

where the unitary matrix [/„ is obtained from IFu and reads 

u 1 ^ -"i) ( 1 

“ detIF„ ; 2(l + u2)-iuo 1 


(6.3) 


(6.4) 


In this indirect way the Kurasov matrices determining the matching conditions that define the 5-5' in¬ 
teractions are related to a subset of the U{2) group which in turn characterizes a variety of self-adjoint 
extensions of the Hq operator in the standard manner. 

A remarkable fact is the following: even though My and IFu are singular matrices at the decoupling 
limit ui = ±1 the corresponding unitary matrices are regular: 


Uy 


i;i = l 


0 -1 

4+ivo n 

A—ivo 



n 4+ivo 

4—ivQ 
-1 0 


(6.5) 
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Writing the linear system (6.4) for the decoupling limit of the coupling vi = 1 we obtain two equations. 
First, 

^ ‘/?( 9 ")= 0 , ( 6 . 6 ) 

i.e., Dirichlet boundary conditions are satisfied reaching the point q from the left. Second, 




4 — itiQ 


((/7(g+) - + i<f'{q^) 


- = 0 

Vq 


(6.7) 


sets Robin boundary conditions at q coming from the right. If the other decoupling value for the coupling, 
vi = —1, is chosen the situation is identical but the conditions at q from the left or from the right are 
exchanged. The results above is in perfect agreement with the Kurasov matching conditions at decoupling 


limit of the couplings, as expressed in formula (5.2). Finally, we mention that in the usual treatment this 


situation is determined from diagonal rather than anti-diagonal matrices. To meet this criterion one merely 


multiply (6.5) by the Ui Pauli matrix. 


Finally, we extract some physical consequences from the non-abelian superposition law for the S-S' 
potentials: 


1. The quantum self-energy of two S-S' configurations due to quantum vacuum scalar fluctuations should 
inherit somehow some features from the non-abelian composition law. A wel known procedure 
to regularize such divergent quantity is to start from the heat trace of the Hamiltonian operator: 
hnit) = Tri 2 exp[—ti7], where t is Schwinger proper time. This spectral function is obtained for any 
Schrodinger Hamiltonian in terms of the bound state energies and the spectral density, defined in 
turn from the total phase shift, see the two formulas just above (32) in [^. From the total phase 
shift produced by the S-S' potential determined by the Kurasov matrix and the bound state 
energy we hnd: 


hff {t- My) = ^2i7rerfc ^ 2 ^ ^ 

where erfc is the complementary error function and 9 the Heaviside step function. It is of note that the 
exact heat kernel for a single S-S' is well defined in the decoupling limit Vi —>■ ±1. We conjecture that 
taking the zero distance limit in the double S-S' potential, the heat trace for the superposed potential 
will be of the form hn (t; My ■ My,) = hn {t; My). The quantum vacuum energy is essentially the 
spectral zeta function evaluated at s = —1/2. This second spectral function is obtained from the 
heat trace via Mellin’s transform of the heat trace: Ch{s) = hnit), where r(s) is 

Euler Gamma function and s a complex parameter, s = —1/2 is a pole of the meromorphic function 
Cnis) in C. One regularizes the divergent quantum vacuum energy by assigning to Ci? its value in 
a regular point. What the non-abelian addition law tell us are the values of the parameters uq and 
Ml as functions of vo,wo,vi, wi entering in this regularized expression after the q = 0 limit has been 
taken. 

2. An ideal model of electric conductivity in solids is provided by a ^-Dirac comb where 5-point interac¬ 
tions sit at the ions sites. This model can be enriched in a twofold way: (1) S' potentials are added 
at every site of the lattice. (2) Two species of ions, henceforth two species of S-S' interactions, such 
that the solid is characterized by the periodic potential 

V(x) = (voS(x — nq) -\- viS{x — nq) -\- uoS{x — nq -\- p) uiS(x — nq-\- p)) (6-9) 


are considered. There are two limits to a single species: p ^ 0,q. Unlike in the Dirac comb of two 
species where the two limits are equivalent, the two merging processes are different in the (6.9) comb 
due to the non-abelian superposition law. In fact, allowing p to vary we have an infinitely repeated 
Casimir piston, which can be reduced to one piston by restricting the system to the primitive cell. 
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